HEAT KERNEL ASYMPTOTICS FOR LAPLACE 
TYPE OPERATORS AND MATRIX KDV 
HIERARCHY 



Abstract. We study the heat kernel asymptotics for the Lapla- 
ce type differential operators on vector bundles over Riemannian 
manifolds. In particular this includes the case of the Laplacians 
acting on differential p-forms. We extend our results obtained 
earlier for the scalar Laplacian and present closed formulas for 
all heat invariants associated with these operators. As another 
application, we present new explicit formulas for the matrix 
Korteweg-de Vries hierarchy. 



1. Introduction and main results 

1.1. Laplace type operators. Let M be a compact smooth d- 
dimensional Riemannian manifold without boundary with a metric 
(gij) and V be a smooth r-dimensional vector bundle over M. Let 
(x±, . . . ,Xd) be a system of local coordinates centered at the point 
x G M and (v i, . . . , v r ) be a local frame for V defined near the ori- 
gin. 

Consider a Laplace type differential operator, i.e. a second-order 
elliptic self-adjoint differential operator D : C°°(V) — > C°°(V) with 
the scalar leading symbol given by the metric tensor (see [Gl]). In 
the chosen local system it can be written as 
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where (g 1 ^) denotes the inverse of the matrix (g^), / is the iden- 
tity matrix and = Bk(x\, . . . ,Xd), C = C(xi, . . . ,Xd) are some 
endomorphisms of the bundle V given by r x r square matrices. 

In particular, any Laplacian A p acting on the ^-vector bundle 
Q P (M) of differential p-forms can be presented in this way (see [Ga]). 
Another example of a Laplace-type operator is a (1-dimensional) ma- 
trix Schrodinger operator (see section 3.1): 

L = I-^ + U(x), (1.1.2) 
where U (x) is a r x r hermitian matrix. 

1.2. Heat kernel asymptotics. Let us define the heat operator 
e -tD f or ^ > o (see [Gl]). It is an infinitely smoothing opera- 
tor from L 2 (V) — > C°°(V) and it is given by the kernel function 
K(x, y, t, D) : V y — > V x as follows: 

(e- w f)(x) = J K(x,y,t,D)f(y)Jgdy. 

M 

where g = det(gij). The function K(x, y, t, D) is called the heat kernel 
associated with the operator D. It has the following asymptotic 
expansion on the diagonal as t — > 0+ (see [Se], [Gl]): 

oo 

K(x,x,t,D) ~J2Mx,D)t n ~i 

where A n are certain endomorphisms of the fiber V x . 

Consider the spectral decomposition of the self-adjoing operator D 
into a complete orthonormal basis of eigenf unctions ipi and eigenval- 
ues Aj, < Ai < A 2 .... Then 

K(x, y, t,D)=Y, e- Xit Vi{x) ® ^(y) 

Denoting a n {x, D) = Tr(A n (x, D)) where Tr is the fiber trace of the 
endomorphism A n (x,D), we get 

oo 

Tr{K{x,x,t,D)) =Y,e- tXi {ipi,ipi){x) ~ £ a n (x, D)t n ~l 

i n=0 

The coefficients a n (x, D) are the heat invariants of the Riemannian 
manifold M associated with the operator D. They are homogeneous 
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polynomials of degree 2n in the derivatives of the Riemannian metric 
{g ij } at the point x ([G2]). 

Computation of heat kernel coefficients is a long-standing prob- 
lem in spectral geometry (we refer to [P3] for a historical review). 
In [P3] we calculated heat invariants of Riemannian manifolds for 
the scalar Laplacian using the Agmon-Kannai (commutator) method 
(see [PI]). In this paper we present its generalization and extend 
our formulas for the heat invariants obtained in [P3] to the general 
Laplace type operators. As another application, we extend the ex- 
plicit expressions for the scalar Korteweg-de Vries hierarchy obtained 
in [P2], [P3] to the matrix case using the heat kernel asymptotics of 
the matrix Schrodinger operator. 

The main difficulty in applying the Agmon-Kannai method to 
matrix-valued operators is that the sum of the orders of two such 
operators is not necessarily greater than the order of their commu- 
tator (cf. fl2.2.1|) ). However, for Laplace type operators this problem 
in fact does not appear (see section 2.2). 

1.3. Main result. Denote by p x : M — > R the distance function on 
the manifold M: for every y e M the distance between the points y 
and x is p x (y)- Let us also consider the operator D and its powers as 
r xr matrices (cf. (|1 . 1 . 1| ) ) . Let Tr(D^ +n ) denote the scalar differential 
operator which is the matrix trace, i.e. the sum of diagonal entries 
of the operator D^ +n . 

Theorem 1.3.1. The endomorphisms A n (x,D) are equal to 



A n (x,D) = (47r)- rf / 2 (-l)" £ )u '/^ (D> +n )( P M 2j • *) 

i=o J - \J > n )- 



3/i 




and hence the heat invariants a n (x,D) are given by 

3n (3n+d/2\ 




y=x 



(3n + d/2\ _ T(3n + 1 + d/2) 



(1.3.2) 



\j + d/2) (3n-j)\T(j + l + d/2Y 



Theorem |1.3.1| is proved in section 3.1. 
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1.4. Remarks. The endomorphisms A n (x, D) do not depend on the 
choice of the coordinate system on M and the local frame for V. 
However, expressions in Theorem |1.3.1| depend on the frame in the 
fiber V x . This means that these expressions are subject to certain 
combinatorial cancellations which should reveal the invariant nature 
of the endomorphisms A n (x,D). Note that independence on the 
choice of the coordinate system is obtained exactly through the study 
of the combinatorial structure of our formulas (see [P3]). 

Another indication of such hidden cancellations (as was noticed by 
P. Gilkey) is the explicit dependence of the expressions in Theorem 
|1.3.1| on the dimension d of the manifold M. Indeed, if one writes 
A n (x,D) in terms of the curvature tensor and its derivatives, the 
dimension never appears. Therefore, there should exist a way to 
rewrite Theorem |1.3.1| in a dimension-independent form as well. 



2. The Agmon-Kannai method for Laplace type 

operators 

2.1. The Agmon-Kannai expansion. The method for computa- 
tion of heat invariants developed in [PI], [P3] is based on an as- 
ymptotic expansion for resolvent kernels of elliptic operators due to 
S. Agmon and Y. Kannai ([AK]). Let us briefly review the Agmon- 
Kannai result. 

Definition 2.1.1. (see [AK]). Let J = (ji, . . . ,ji) be a finite vector 
with nonnegative integer components and let P, Q be linear operators 
on a linear space M. The multiple commutator [P, Q; J] is defined 
for all such vectors J by 

[P,Q;J] = (adP)*Q---(adP)*Q, 
where (ad P)Q = [P,Q\. 

Let if be a a self-adjoint elliptic differential operator of order s on 
a manifold M of dimension d < s, and Hq be the operator obtained 
by freezing the coefficients of the principal part H' of the operator 
H at some point x G M: H = H'(x). Denote by R\(x, y) the kernel 
of the resolvent R\ = (H — A) -1 , and by F\(x,y) — the kernel of 
Fx = (Ho — A) -1 . The Agmon-Kannai formula in its original form 
reads as follows ([AK]): 
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Theorem 2.1.2. The following asymptotic representation on the di- 
agonal holds for the kernel of the resolvent R\ as X — ► oo: 



1 oo 

Rx(x,x) ~ — (F x (x,x) + £([Pf -H,H ;J}F^ +r+1 (x,x))), 



V9 



(2.1.3) 



where the sum is taken over all vectors J of length > 1 with nonneg- 
ative integer entries. 

2.2. An extension of the Agmon-Kannai theorem. Denote by 
o(P) the order of the differential operator P (if P = (Py) is a matrix 
operator then o(P) = maxy o(Py)). For any two scalar differential 
operators P, Q we have: 

o([P,Q]) <o(P) + o(Q)-l (2.2.1) 

This simple inequality lies in the basis of the Agmon-Kannai commu- 
tator expansions (cf. [AK], Lemma 5.1). Clearly, in the matrix case 
the relation ( 2.2.1|) does not hold in general, and therefore Theorem 
2.1.2| is no longer true. However, under certain assumptions Theorem 
[2.1. 2| can be extended to the matrix case in a straightforward way 
(see [Kan] for an analogue of this theorem in a much more general 
situation) . 

Lemma 2.2.2. Let H be anrxr matrix differential operator of order 
s with a scalar principal part H' . Then the asymptotic expansion 
( ft.l.yp remains valid. 

Proof. Indeed, since H' is a scalar operator, Hq = ho ■ I is also 
scalar, o(Hq) = s. Let P be an r x r matrix differential of order p. 
Consider the commutator [P, Hq], Since H is scalar, the commutator 
is a matrix operator with the entries [Py, ho]. For each entry we have: 
o([P i:j , h ]) < p+s-1, i,j = 1, . . . ,r. Therefore, o ([P,H ]) < p+s-1. 
Let us note, that all commutators in the formula ( |2 . 1 . 3| ) are exactly of 
the form [P, Hq] , where P is some matrix operator. We have shown 
that inequality ( |2.2.1| ) holds for all commutators of this kind, and 
therefore in this case we may just repeat the proof of Theorem |2.1.2 
(see [AK]). This completes the proof of the lemma. □ 

In [PI] we found a concise reformulation of the Agmon-Kannai 
theorem for scalar operators (Theorem 1.2). Lemma |2.2.2| allows to 
extend it immediately to the case of matrix operators with a scalar 
principal part. 
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Theorem 2.2.3. Let H be an elliptic differential operator of order s 
acting on a r- dimensional vector bundle V over a Riemannian man- 
ifold M of order d < s. The following asymptotic representation 
on the diagonal holds for the resolvent kernel R\(x,y) of an elliptic 
operator H as A — ► oo : 

1 oo 

Rx(x,x) ~ — J2 X m F™ +1 (x,x), (2.2.4) 
V 9 m =o 

where the operators X m are defined by: 

Note that the condition d < s can be relaxed in the same way as 
in [P3] - - by taking the derivatives of the resolvent with respect to 
the spectral parameter A (see Theorem 2.3.1 in [P3]). 

2.3. Proof of Theorem 1.3.1. Indeed, consider a normal coordi- 
nate system centered at the point x G M. Then g%j{x) = Sij and 
hence the principal part of the operator D with the coefficients frozen 
at the point x has the form D = A • I, where A is the scalar Lapla- 
cian on R d and I is the r x r identity matrix. Due to Theorem 2.2.3| 



we can repeat the arguments of the proofs of Theorems 1.2.1 and 
3.1.1 from [P3] just taking into account multiplication by the matrix 
I. Finally we obtain the expression for A n (x, D) identical to the for- 
mula (4.2.2) in [P3]. Expressing the distance function in terms of 
normal coordinates as in [P3] we obtain the formulas for the endo- 
morphisms A n (x, D). Taking the traces of A n (x, D) we get the heat 
invariants a n (x, D) and this completes the proof of the theorem. □ 



3. Matrix KdV hierarchy 

3.1. Heat kernel asymptotics for Schrodinger operator. Con- 
sider the 1-dimensional matrix Schrodinger operator L (see (|1.1.2| )) 
with a potential U which is an hermitian r x r matrix. Its heat kernel 
H(t,x,y) is the fundamental solution of the heat equation 
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It has the following asymptotic representation on the diagonal as 
t -> 0+: 

1 oo 

H(t,x,x)~—=JTh n [U]t n , 

where /i n [C/] are some polynomials in the matrix U(x) and its deriva- 
tives. 

The matrix KdV hierarchy is defined by (see [AvScl]): 

dU d 

m = »?>M (3L1) 



where 



G n [U] = ^h n [U], neN. 



Set Uq = U, U n = d n U/dx n , n G N, where U n , n > are formal 
variables. The sequence of polynomials G n [U] = G n [U , U\, U 2 , ■ ■ ■] 
starts with (see [AvSc2]): 

Gi[U) = U Q , G 2 [U]=U 2 + 3Ul 

G 3 [u] =U 4 + 5U U 2 + 5U 2 U + Wl + 10f/ 3 , . . . 

3.2. Computation of matrix KdV hierarchy. In [P2], [P3] we 
presented explicit formulas for the scalar Korteweg-de Vries hierarchy 
using heat invariants of the scalar 1-dimensional Schrodinger operator 
(we refer to [P2] for the history of this question). Theorem |1.3.1| 
allows to extend our results to the matrix KdV hierarchy as well. 

For other approaches to explicit computations of the matrix KdV 
hierarchy see [AvScl], [AvSc2]. 

Theorem 3.2.1. The KdV hierarchy is given by: 

where the polynomial P n j[U] is obtained from L :j+n (x 2: ')\ x=0 by a for- 
mal change of variables: £/j(0) — > U i} i = 0, .., 2n + 2j — 2. 

This expression can be completely expanded due to a formula for 
the powers of the Schrodinger operator ([Rid]) which remains valid 
in the matrix well (note that now Ui are matrices and hence 

do not commute). 
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Theorem 3.2.2. The polynomials G n [U], n G N are equal to: 

(2,,.)! ^(n+j\ tmm ^ 

fc^-j hfcp — 2{n — p) 

c k v - M (n p -i - h - ■ ■ 

0<lo<h<-<l P -i=j+n-p ^2 / V 

2Zi>fciH hfci+i, i=0,...,p-l. 
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